arXiv:1501.03014vl [astro-ph.CO] 13 Jan 2015 


Astronomy & Astrophysics manuscript no. Camelio_Lombardi ©ESO 2015 

January 14, 2015 


On the origin of intrinsic alignment in cosmic shear 
measurements: an analytic argument 

Giovanni Camel i(P-d and Marco LombardPS 

1 University of Milan, Department of Physics, via Celoria 16,1-20133 Milan, Italy 

2 University of Rome “La Sapienza”, Department of Physics, p.le Aldo Moro 2,1-00185 Rome, Italy 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

Received ***date***; Accepted ***date*** 


ABSTRACT 

Galaxy intrinsic alignment can be a severe source of error in weak-lensing studies. The problem has been widely studied by numerical 
simulations and with heuristic models, but without a clear theoretical justification of its origin and amplitude. In particular, it is still 
unclear whether intrinsic alignment of galaxies is dominated by formation and accretion processes or by the effects of the instantaneous 
tidal field acting upon them. We investigate this question by developing a simple model of intrinsic alignment for elliptical galaxies, 
based on the instantaneous tidal field. Making use of the galaxy stellar distribution function, we estimate the intrinsic alignment signal 
and find that although it has the expected dependence on the tidal field, it is too weak to account for the observed signal. This is an 
indirect validation of the standard view that intrinsic alignment is caused by formation and/or accretion processes. 
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1. Introduction 

It has long been recognized that gravitational lensing is a partic¬ 
ularly robust method to investigate the mass distribution of mas¬ 
sive astronomical o bjects (e.g., see Refsdal 1 964tlFu et al.l l2008: 
Jauza c et aj]l2012h . Gravitational lensing is insensitive to the 
chemical and physical state of the deflecting matter and there¬ 
fore treats ordinary and dark matter components of the deflector 
squally. In addition, it depends on a clearly understood physi¬ 
cal process, the distortion of the space-time induced by masses 
as described by General Relativity, for which it is not required 
that the deflector be in a state of dynamical equilibrium. As a 
result, it currently represents our best tool to study the matter 
distribution in the Universe. 

A particularly important use of weak-lensing methods, 
which are based on the statistical analysis of the weak corre¬ 
lation among the observed shapes of distant galaxies, is studying 
the large-scale structure of the Universe. In weak cosmologi¬ 
cal lensing, the deflecting matter studied (in contrast to the case 
of a cluster of galaxies) is distributed along the entire line of 
sight. As a result, the distinction between source and lens is less 
clear: the light of every galaxy that contributes to the distortion 
of space-time is bent by the distribution of matter along its path. 
The correlation between the apparent ellipticities of every pair 
of galaxies determines the signal; this is called cosmic shear. 
Since cosmic shear is related to the power spectrum of th e den¬ 
sity contrast (e.g., see Miral da-Escudel 199 It iKaiserll 1 9921) . with 
this tool one measures the statistical properties of the large-scale 
structure in the Universe. Therefore, weak cosmological lensing 
is an independent method for testing cosmological models (e.g., 
see iKilbinger et all20l3h . 

Cosmic-shear measurements are particularly challenging: 
the signal is very weak and is affected by many error sources. 
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From the modeling point of view, one of the main difficulties 
is related to the so-called intrinsic alignment, that is, the galaxy 
intrinsic shape alignment, which is not caused by gravitational 
lensing, but by the source gravitati onal field ( Heave ns et al.l 
2000: (Croft & Metzlerll200(il : see also lTroxel & Ishakll2014l for~a 
recent review). Although the effects of intrinsic alignment can be 
reduced by removing pairs of galaxies that are physically close 
to the measurement of the shear correlation dKing & Schneiderl 
l2002h IHevmans &~H eavens 2003), there remains a residual bias 
(called GI mode, see below). Therefore, it is important to clar¬ 
ify these effects and to assess their impact on current and fu¬ 
ture cosmic-shear surveys. Although some heuristic models 
for intrinsic alignmen t are available (e.g.. ICatelan et all 120011 : 
iHirata & Seliakl 120041) . they are unfortunately incomplete; in 
particular, they lack an analytic estimate of the bias magnitude 
(in general, such an estimate would require the use of compli¬ 
cated galaxy formation and stellar dynamical models). There¬ 
fore, one has to resort to the bia s magnitude measured from sim¬ 
ulations dCroft & Metzledl2000t He avens et aT]l2000tliing||2002l : 


IHevmans et alJl200flt Tenneti et al .1120141) for a comparison with 

the observations. 


The exact origin of intrinsic alignment is still debated (e.g., 
see lPereira & Kuhnl2005h . and it is still unclear whether intrinsic 
alignment of galaxies is dominated by formation and accretion 
processes or by the effects of the instantaneous tidal field acting 
upon them. The commonly accepted view, supported by simu¬ 
lations and observations, is that the intrinsic alignment is ca used 
by the tidal field acting at the formation epoch (Catela n et al.l 
1200ll:fHirata & Seliakil20041) and possibly on the merger history. 
In this paper we strengthen this view with analytic arguments 
using a reductio ad absurdum: that is, we start with the assump¬ 
tion that intrinsic alignment is dominated by the instantaneous 
tidal field currently at work on galaxies and show that the pre¬ 
dicted effect is several orders of magnitude weaker than the ob- 
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served effect. In our model, we consider a spherical early-type 
galaxy and study the deformation of its shape to first order after 
introducing a weak external tidal field. We describe the effects 
of the external tidal field on the intrinsic (i.e., unlensed) shape 
of the elliptical galaxy, determining the luminous quadrupole 
of the galaxy by means of its collisionless stellar distribution 
function. The galaxy ellipticity is then computed through the 
luminous quadrupole. We find the expected dependence of the 
ellipticity to the tidal field, as postulated by the heuristic model 
of ldatelan et al . ( 200 11. but the constant of proportionality deter¬ 
mined by our model is more than four orders of magnitude lower 
than that observed dJoachimi et al.ll20TTh . and consequently can¬ 
not account for the intrinsic alignment. 

The paper is organized as follows: in Sect. [2] we present the 
problem of intrinsic alignment, summarize the relevant classi¬ 
fications, and outline the properties of different types of align¬ 
ment. In Sect. 0 we introduce the stellar distribution function 
and estimate the ellipticity induced by an external tidal field on 
an elliptical galaxy. In Sect. 0 we compare our result with that 
of lJoachimi et ali (1201 ll) and apply our model to the galaxy dark 
matter halo. We draw our conclusions in Sect. [5] In Appendix lAl 
we study the dynamics of elliptical and spiral galaxies in the 
rigid-body approximation. Finally, in Appendix[B] we derive an 
alternative intrinsic alignment ellipticity determined by instanta¬ 
neous tidal fields by means of the equipotential surfaces of ellip¬ 
tical galaxies. Throughout this paper we take a standard cosmo¬ 
logical model with Hubble constant Hq = TOkms' 1 kpc~ L , mass 
density parameter Q m =0.3. 


2. General context 

2.1. Cosmic shear 


In principle, unbiased shear measurements can be obtained from 
the luminous quadrupole tensor Q, , of each galaxy in a given 
field of view (e.g JBartelmann & Schneideifl200ll) 


f I(e)didj d 2 d 
'" J ~ f 1 ( 6 ) d 2 8 


( 1 ) 


where 1 ( 6 ) is the surface brightness observed at angular posi¬ 
tion 6 , defined as the two-dimensional vector from the luminous 
center of the galaxy to a given point. In reality, this definition 
of quadrupole moments produces measurements with vanishing 
signal-to-noise ratio because the noise at large radii 6 = ||0|| dom¬ 
inates Eq. (0. For this reason, in actual weak-lensing observa¬ 
tions one resorts to weighted quadrupole tensors or to alternative 
measurement techniques (e.g., iKaiser et alj|1995l) . 

For simplicity, we use definition (□) here, which 
has convenient transformation properties. In this case 
(ISchramm & KavseHl 19951 ISeitz & SchneideJl 19971) and under 
the assumption that the unlensed source ellipticities have van¬ 
ishing average, the observed complex ellipticity, defined as 


where am and a m are the semi-major and semi-minor axes of 
the isophotes. If the light paths are bent weakly, the observed 
ellipticity s of the galaxy is (e.g. jBartelmann & Schneiderl200li) 

£ = e s + y , (4) 

where s s is the intrinsic (i.e., unlensed) shape of the galaxy and 
y is the cosmic shear. 

In cosmic-shear studies it is convenient to define three ellip¬ 
ticity correlation functions: 

C ++ (6) = <£ + (0 o )£ + (0o + 0))e o , (5) 

Cxx(6) = (Sx(6o)s x (0o + 6))g 0 , (6) 

C +x (6 ) = (s + (0 o )sx(0o + 0))e o , (7) 

where s+ and s x are the real and imaginary part of the complex 

ellipticity in a coordinate system with the real axis along the 
line connecting the two correlating points (the centers of the two 
galaxies, 6q and 6o + 6). The correlation functions only depend 
on the modulus of 6 because of isotropy. 

We now consider a pair of galaxies and use the superscripts 
fg and bg to denote foreground and background objects. Then 
we find 


Cij = (sfs b h = (yfy^+y[ 


fg s,bg 


s.fg bg , s,fg s,bg\ 

■ S "' S + £, S £. S ) 


+ s i°y 
<7,GG 


' 1 I J 

= Cii.GG + Gy, GI + G 


d,n 


( 8 ) 


where [i; j] = {+; x), y fg (y bg ) is the cosmic shear acting on the 
foreground (background) object, C,jx;,G = (y}y b f) is the cos¬ 
mic shear correlation term, and G/gi = 0'[ s £*’ bs + £;' fs y bs ) and 

Gy ,11 = (£-’ fs £y’ bs ) are the gravitational-intrinsic and the intrinsic- 
intrinsic correlation terms (see below). We note that the dom¬ 
inant term in Gy.oi is ( e / fS 7y S )i the other term, (y[ S £*’ bg ), is 
only expected to be non-negligible if there are long-range cor¬ 
relations in the gravitational tidal field (e.g., a dark matter fila¬ 
ment along the line of sight), and can be ignored in most situa¬ 
tions. In the case of a flat cosmology, the correlation functions 
are related to the density contrast power spectrum P$ (e.g., see 
iBartelmann & SchneideH[200lllCatelan et al.ll2001h 


C++, Gg(0) + G< x ,gg(#) - 

f°° Idl 

J ’ 

(9) 

C++, Gc(0) - Cxx,G g(#) = 

r°° (d{ 

J 4(m ’ 

(10) 

C+x,Gg(#) = 

0 , 

(11) 

Py(t) = 

9 H*ni pH W 2 (w ) | 


4c 4 Jo a 2 (w) °' 



(12) 

W(w) = 

f“’ H , w' -w , 

P (w ) - - — dw . 

Jw W 

(13) 


Qn - Q22 + 2zQi2 

£ = - ,- , ( 2 ) 

Tr Q + 2 VdetQ 

provides an unbiased estimate of the (reduced) shear. For a lumi¬ 
nosity distribution with constant flattening elliptical isophotes, 
the modulus of the complex ellipticity is 

■ | a m — 

M =-, (3) 

Am + Am 
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Here P y is the cosmic-shear power spectrum, Jq and J 4 are 
Bessel functions of the first kind, p(w) is the probability of de¬ 
tecting a galaxy at comoving distance w, and u;h is the comoving 
horizon distance. 

To obtain relevant cosmological results from cosmic-shear 
measurements, it is important to understand each error source. 
From the theoretical point of view, an important error source 
is the contribution to the ellipticity correlation induced by the 
intrinsic alignment, Gy.Gi + Gy.n- 
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2.2. Intrinsic alignment 
As noted previously in computational dHeavens et al. 2000; 


Croft & Metzled l2C)00t) and theoretical studies ( 1C ate I an et al.l 
200 lib unlensed ellinticities need not be distributed isotropically. 
In particular, close pairs of galaxies may exhibit correlated ellip- 
ticities as a result of their mutual gravitational interaction and of 
an external gravitational potential. 

The intrinsic alignment (IA) is the contribution to the ob¬ 
served galaxy ellipticity correlation, caused by the gravitational 
tidal field in wh ich galaxies are placed. In the literature it is 
usually assumed dCatelan et al.l200iT ; lHirata & Seliakl2004l) that 
the galaxy intrinsic ellipticity is frozen at the formation epoch. 
IA is_dassified Recording to the physical generating mechanism 
(ICatelan et alil20071) . which is different according to the type of 
galaxy that the external gravitational field acts upon. 

In elliptical gala xies, the external gravitational field stretches 
the galaxy shapes (ICroft & Metzleii 120001) . partly determining 
their intrinsic ellipticitiesQ In contrast, in spiral galaxies the in¬ 
trinsic ellipticity is related to their angular momentum, produced 
by the torque provided by the external gravitational field during 
galaxy formation and by projection effects dHeavens et al .1200(1) . 

From a geometrical point of view dHirata & Selia k 2004), 
the intrinsic correlation (for a second-order statistics, like the 
shear two-point correlation function) is made of two contribu¬ 
tions: intrinsic-intrinsic (II) and gravitational-intrinsic (GI). The 
II contribution is the correlation of the unlensed ellipticities of 
two physically clos43 source galaxies, generated by a corre¬ 
lation in the gravitational field felt by the galaxies (which is 
caused by their mutual interaction and by the external large- 
scale structure). In contrast, the GI correlation is the correlation 
between the observed ellipticities of two galaxies that are dis¬ 
tant from each other, but angularly close in the sky view plane 
dHirata & Seliak82004l) . This correlation is caused by foreground 
mass inhomogeneities that produce two effects through the tidal 
field: they induce an intrinsic ellipticity in foreground galaxies 
through direct gravitational interaction and modify the observed 
ellipticities of background galaxies by means of gravitational 
lensing. 


2.2.1. Elliptical galaxies 

Elliptical galaxies are expected to be polarized by the ex¬ 
ternal tidal field they are located in (ICroft & Metzlerl [2000; 
ICatelan et al.ll200TI) . 

4 = C(d] - d\)U cn , ( 14 ) 

4 = 2Cd t d 2 u cn . ( 15 ) 


1 Elliptical galaxies may be non-spherical on their own, independently 
of any tidal stretching because of anisotropic pressure and angular mo¬ 
mentum. 

2 In cosmic-shear measurements, the II signal is usually suppressed 
by excluding pairs of galaxies from the analysis that are thought to be 
physically close (with their physical distance estimated from photomet¬ 
ric redshifts and angular distance in the sky). However, the assumption 
that only physically close galaxies have II alignment might not hold 
because the II alignment could arise also in far away galaxies if there 
are long-range correlations in the gravitational field (e.g., those asso¬ 
ciated with a dark matter filament, the typical length scale of which is 
numerically estimated a s IO()/z 1 Mpc by Srringel et al. 2005; see also 
iMandelbaum et al,l2006l : lLed2004l : lHirata et alj2007l) . Moreover, since 
the distance between galaxies is estimated from photometric redshifts, 
galaxies thought to be physically distant may actually be physically 
close as a result of errors in the photometric redshift estimates. 


Here e s is the galaxy ellipticity in the plane perpendicular to 
the line of sight, f/ ext represents the large-scale potential (at 
the galaxy formation epoch) associated with inhomogeneities 
of cosmological origin, and C is a constant that in principle 
could be determine d by “a complete gala ctosynthesis model” 
(ICatelan et al.ll20071) . ICatelan et al.l (1200 ll) indirectly estimated 
this constant in a heuristic way by means of the relation 


<4) = -C 2 



k 2 P s {k) 


3W d , 

kR I 


( 16 ) 


where R = I h 1 2 Mpc is the characteristic scale over which 
a galaxy forms. An II alignment for elliptical galax¬ 
ies has indeed been detected in numerical simulations 


2005; Agustsson & Brainerd 2006 

H rata et al. 2007 

Faltenbacher et al. 2009; Joachimi et al. 

2011). 


For early-type galaxy samples with broad redshift distribu¬ 
tion and for typical foreground mass inhomogeneities, the GI 
contribution is non-vanishing and is expected to be greater than 
that of the II contribution dHirata & Seliakll2004t) . The GI term 
is due to the tidal effects of a mass overdensity on foreground 
galaxies, and the lensing effects of the same overdensity on back¬ 
ground galaxies. Therefore, the relevant factor is the probability 
of observing an elliptical galaxy integrated along the line of sight 
from the mass overdensity to the horizon. In contrast, the rele¬ 
vant factor for the II term is the probability of observing a second 
galaxy at the same redshift as the first galaxy. 

Finally, we note that according to this model, the GI align¬ 
ment produces an anticorrelation of ellipticities (Cgi < 0): 
the foreground elliptical galaxy physical (intrinsic) shape is 
stretched along the gravitational tidal field, whereas the back¬ 
ground galaxy apparent sh ape is s tretched perpendicular to it 
because of lensing dHirata & Selia k 2004). For example, in a 
galaxy cluster, cluster members would be preferentially aligned 
radially, whereas the background galaxies would be seen prefer¬ 
entially tangentially. 


2.2.2. Spiral galaxies 

In spiral galaxies the II contribution is thought to be caused by 
the torque provided by the external tidal field during formation 

dHeavens et alJl20()7t ICatelan et alJl200lfc iHirata & Seliakll2004l 

12010l) . This contribution is to second order in the external tidal 
field because the tidal field first has to generate an anisotropic 
moment of inertia l, ; (as in elliptical galaxies, the distribution of 
mass is expected to be elonga ted along the gravitational field gra¬ 
dient), and then to torque it dPeebles|] l969t lDoroshkevichll 19701 : 
Whitel 179841) . The torque provided by the external tidal field 
generates the angular momentum of the spiral proto-galaxy (this 
is similar to what we show in Appendix [A] where we compute 
the torque generated by the external tidal field on an elliptical 
galaxy not aligned with it). One then needs to specify the corre¬ 
lation between the angular momentum and the galaxy ellipticity, 
for example by assuming zero thickness for the (circular) disk 
of the galaxy and considering the projection effects, as reported 
bv lHeavens et al.l d2000l) . Finally, since a correlation in the tidal 
field induces a correlation in the angular momenta of close spiral 
galaxy pairs, the expecte d II power spect rum can be computed 
analytically, as shown bv lHirata & Seliak! (120041) . 

The II contribution for spiral galaxies is of higher order 
in the tidal field and thus should be lower than for ellipticals 
dHirata & Seliak|[2004l) . In addition, little GI contribution is ex¬ 
pected for spiral galaxies dHirata & S eliak 2004). At present, 
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there is no clear observational evidence of IA in spiral galaxie s 

dHirata et al.l2007llFaltenbacher et al.l2009tlBlazek et al.ll201 2l) . 

In Appendix lAl we determine the rigid-body precession pe¬ 
riod of spiral galaxies, which is longer than the galaxy deforma¬ 
tion time. This shows that no precession-driven mechanism for 
an alignment of spiral galaxies with the external tidal field can 
be devised. 


3. Distribution function method 

Even if galaxies are many-body systems, it might naively be 
thought that their alignment is driven by rigid-body motions 
(e.g., oscillations and precessions) under the effects of exter¬ 
nal gravitational forces. However, an analysis of the relevant 
time scales shows that rigid-body motions are much slower than 
those ch arac terizing internal dynamics (Appendix [A] see also 
iCiotti & Duttal119941) . This means that at least for IA purposes, 
the galaxy is best studied as a deforming body. In this section, 
we start from the hypothesis that galaxies are not subject to any 
intrinsic alignment during formation and that intrinsic alignment 
is entirely due to tidal fields at the observation time, to which 
galaxies react immediately. As shown in Sect. [4] this assump¬ 
tion leads to an extremely low IA, completely inconsistent with 
observations. This is an indirect proof that IA is driven by for¬ 
mation and merging. 

The deformation of an elliptical galaxy subject to an exter¬ 
nal gravitational field can be modeled in a simplified way by 
studying the deformation of its equipotential surfaces (see Ap¬ 
pendix |B}. and in a more complete manner, by means of the 
stellar distribution function. In both cases, we start with unper¬ 
turbed spherical galaxies for simplicity and study the ellipticity 
induced by the external tidal field. The spherical assumption is 
quite strong, and it may not be straightforward to generalize our 
results to the case of an ensemble average of elliptical early-type 
galaxies. Nevertheless, since we are interested in an order-of- 
magnitude estimate, the spherical assumption is sufficient for our 
purposes. 


3.1. General case 

A galaxy is a complex many-body system, w ith a gravity-driven 
dynamics (IChandrasekhail[l942t lBertinll2000l) . In this context, 
an important tool is the stellar distribution function /*(x, x, t), 
which is, essentially, the stellar density distribution in phase 
space. We also define (e.g.. iBertinl 12000 ) the (total) distribu¬ 
tion function / = /* + /dm as the sum of the stellar distribution 
function and the dark matter distribution function. We assume 
the collisionless Boltzmann equation and the Poisson equation: 


df df . df dU 
dt dxi ' dx{ dxj 


(17) 


V~t/(x, t) = 4nGp(x, t ) = 47rGm 


■I 


//(X.X. 

f' 


t) d 3 i 


47tG/h I /*(x,x, t) d J x + AnGm /dm(x,x, f)<rx , (18) 


where m is the mean mass of a galaxy star (and of a DM “par¬ 
ticle” as well). In the following, we use the King model (iKingi 


1 19661 ) for the stellar distribution functiorQ 


/*(x, x, t) = 



-t/(x)/o-J—||x|| 2 /(2o/) _ p-EtJoj 


»] E < E lr , 
E > E lr , 


(19) 


where E = (7(x) + ^ is the single-star energy and <r, is the 
stellar velocity dispersion, and we start from a spherical unper¬ 
turbed (early-type) galaxy with potential Uq. In other words, in 
the unperturbed galaxy the DM supplies the exact contribution 
to have the given potential Uq and the King distribution function 
for/*. 

We take an external tidal potential and add it to the (unper¬ 
turbed) galaxy potential, thus ignoring the changes in the galaxy 
potential induced by the deformation of the galaxy (for a similar 
not self-consistent a pproach, see ICiotti & Duttal 1 1 994t s ee als o 
iBertin & Varrill2008l) . As pointed out bv ICiotti & Duttal (1 19941) . 
it is possible to neglect the actual self-force of the galaxy be¬ 
cause it decreases from the center of the galaxy (see, e.g., the 
potential (l28l > of the unperturbed self-field), whereas the exter- 
nal tidal force increases. Therefore, we take U - Uo + (Aidai (see 
IBertin & Varril 120081: IVarri & Bertinl 120091 for a self-consistent 
approach applied to a globular stellar cluster). We take a weak 
external tidal field |l/ t idai/o"?l ^ 1 and expand the stellar distri¬ 
bution function with respect to it 

/* = - ofOtjXiXj + 

= fT + f* ) + • ■ • > 

y<0) _ A ^t/ 0 (x)/ojHl*ll 2 /(2/b _ e -£,, 

/< 1} - - afOijXiXjA^^-^^) t 

where <1>,y = \djdjU ext is the tidal tensor and U t ;d a i = 
\xjXjdjdjU ext = Q>ijXiXj is the tidal field (potential). 

This expansion is valid only inside the galaxy 
(IBertin & Varril 12008 0: in fact, even if the exponential can 
always be expanded (when |t/ t idai/0"ul <s 1), this expansion is 
not granted to be meaningful unless the perturbed self-potential 
Go + U tidal is not too different from the actual galaxy self¬ 
potential, U. This condition breaks down near the galactic 
boundaries because, as stated above, the galaxy self-potential 
decreases toward the boundaries, whereas the external tidal field 
increases (since Cis traceless, the isopotential surfaces may 
not even be ellipses when | C7tidal I — |Go|). From Eq. <|T]), chang¬ 
ing integration variables to those of standard two-dimensional 
vectors, considering that 7(x) oc j p*(x) djn, and using Eq. ( 1 1 iSl i. 
we obtain for the luminous quadrupole 


...)-Ae- £tt/ ^ 

( 20 ) 

( 21 ) 

( 22 ) 


= Qf) + Q<9 + ... , (23) 

where D s is the angular-diameter distance to the galaxy and N 
is the number of stars of the galaxy. The integration of Eq. (l23l > 
has to be carried out (i) in the region 

Pll ^ Dmax(r) = yj2(E a - U 0 (r )) (24) 

3 The King model is usually applied to study globular clusters; more 
complex distribution functions are generally used for galaxies. Here, 
we choose the King model for its simplicity and ability to clarify the 
key aspects of the dynamical model we are interested in. 
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because of the energy truncation in the distribution function 
and fii) for ||x|| < r max < r tr , a condition that mimics a real 
weak-lensing measurement, where the integration in the lumi¬ 
nous quadrupole is carried out based on a window function. The 
quantity r max is chosen so as to avoid regions where the exter¬ 
nal tidal field becomes similar to or larger than the galactic field, 
that is, the condition |f/tidai(Jmax)l |(/o(r , max)l holds (see discus¬ 
sion above). The zeroth order of the expansion of the luminous 
quadrupole Qjy* is proportional to the identity matrix because the 
unperturbed galaxy field is spherical: 


pf max . 

T\ = I d rr*e- Uolr)ltr ’ I 

Jo Jo 


max _ ) _ _ 

A p -U 0 (r)lo* I y 2 e -a 2 /( 2 ^) d[) 


Jo 


r 4 Q -Uo(r)/(rf 


- Vm*x(r)e t-WW 


\n ,/t>ma x(r) 

+ <W- erf —p— 
2 l V2 tr B 


T 2 = 


-Eu/tr] 


X r max 

r 4 (t»maxW) 


dr, 

3 

dr . 


(25) 


(26) 

(27) 



Fig. 1. Self-potential ( 128b and total mass distribution ( 129b of the 
unperturbed galaxy. The values used are r core = 1 kpc, r tr = 25kpc and 
M = 10 u Ms. The dashed vertical line indicates the truncation radius 
rtr- 


From Eq. (0, considering that Q 4I) oc 5ij (spherical galaxy), we 
find for the complex ellipticity 


We use a galaxy logarithmic potential f/o (Fig. QJ to model the 
total mass distribution in a galaxy (lumino us plus dark matter; 
see e.g. lBertinl2000'[lKoonmans et al.ll2006li 


U 0 (x) = 


po(x) = 


MG 


7*tr ^core 


-MG 


log 

1 

rtr 


/ rcore "b ||x|| 
\ F core 4" ^tr 


Mr tr + ^core core + | |x| | 

' 471 4 l|x||(r CO re + I|X||) 2 

0 


if l|x|| < r,r , 
if l|x|| > r tt , 

if l|x|| < r tr , 
if l|x|| > ^ , 


(28) 


(29) 


where M = 10 n M o is the galaxy mass, Uo(r a ) = 0, r core 
1 kpc is the core radius, and r tr = 25 kpc is the truncation ra¬ 
dius. Throughout this paper we take for an elliptical galaxy 
<x„ =; 200 km s' 1 and r max = 10 kpc (see discussion above). We 
then obtain 


T\ = 6.1 x 10” 5 Mpc 5 (km/s) 3 , (30) 

Ti - 4.6 x 10” 5 Mpc 5 (km/s) 3 . (31) 


The first-order perturbation to the luminous quadrupole is 



A®/,,, r 

NDlal J 

(4 jifATi 
15 crlNDl 


y i y r r ,p-t'o(x)/of-||x|p/(2o^) r |3 
■ ( 2 a>, 7 + f), ; I’r , 


d 3 i 


(32) 


where Trd> = On + O 22 + $ 33 , and 


_ 8 .. Qn-Q^ + ^Qn 

Tr Q (0) + 2 Vdet Q (0) 

Finally, we obtain the general expression for the ellipticity in¬ 
duced on a particular galaxy by an external tidal field. This ex¬ 
pression can be factorized in two parts: one part only depends 
on the particular galaxy we are considering (its mass, its size, its 
self-field, and the velocity dispersion of its stars); the other part 
only depends on the tidal field, with the expecte d dependence 
(e.g., Caitelan et all 1200 lh . The difference with ICatelan et al.1 
(2001) is that in Eq. (l34l ) the external field is that acting at the 
moment of the light emission, and not that acting at the forma¬ 
tion epoch. Comparing Eq. ( l34l) with Eqs. ( 1 1 4b and (fl5l >. we 
obtain the analytic expression for the constant of proportionality 
appearing in Eqs. (fT~4l) and (031 ) 


C = — 


n 


20cr2 Ti - Ti 


=* -4.2 x 10 l4 yr 2 


(35) 


Equation (l34l > is the main result of this paper: it provides a direct 
way to estimate the intrinsic ellipticity of an early-type galaxy 
that is subject to an external tidal field. 

To better appreciate this result, consider a thin lens such as a 
galaxy cluster. In this case, we can write the lens shear as 


7 = c 2 ^^-{d u - d 22 + 2idn) j ~dx 3 t/(x) 

- c~ 2 ls 1 cl ( < f ) ii - <1>22 + 2;d>i2) , (36) 


n = oj r 

Jo 


r 6 Q -U 0 (r)/crl 




/2T Umax(r) 

+ ct„yI- erf 


V2cr„ 

4.9 x 10 -9 Mpc 7 (km/s) 3 . 


dr 


(33) 


where the derivatives are given with respect to the physical co¬ 
ordinates, D s , D\ and Z>i s are the angular distance of a source 
galaxy, of the lens, and between the lens and the source galaxy, 
respectively, and L Q \ is the galaxy cluster typical size. Interest¬ 
ingly, the expression for y has the same form as Eq. (l34l >. that 
is, it depends on the same combination of second-order partial 
derivatives of the tidal tensor <I>. The constants involved, instead. 
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Fig. 2. Radial dependence of the intrinsic galaxy ellipticity in the 
particular cases (Keplerian and DM filament fields) considered in this 
paper, computed through the equipotential approximation (equi, see 
Appendix G3 and through the distribution function formalism (f, see 
Sect. [3}. Note the different radial dependence of the Keplerian (R~ 3 ) 
and the DM filament (R~ 2 ) cases; the equipotential approximation and 
the distribution function formalism yield the same radial dependence in 
the two cases, but different multiplicative factors (see the text for de¬ 
tails). 


where R is the distance of the galaxy from the DM straight line 
(the filament), 6 = diag( 1,1,0) and A is the linear mass dis¬ 
tribution of the dark matter filament. Obviously, aa.tidau = 0. 
We now need to estimate the quantity A , the linear density of 
the DM filament. DM filaments are thought to extend between 
massive galaxy clusters (e.g.. Snrinsiel et alJl2005l) . We there¬ 
fore consider the number density of massive clusters and esti¬ 
mate the average distance d c \ between two clusters. To do so, we 
use the cumulative distribution of galaxy clusters measured by 
iBahcall &Cer3 (119931) : 

n(> M c i) = 4 x 10 5 -e "* h 3 Mpc 3 , (42) 

M ci 

d c i(M cl ) = (n{> M cl ))~ 1/3 , (43) 

where M* = (1.8 ± 0.3) x 10 14 /i 1 M Q and h = 0.7. Setting 
M cl = M* leadfQto 

d c \ =* 58 Mpc , (44) 

which is consistent with the typical value reported by 
ISnrinuel et al.l (120051) . that is, 100Mpc. We assume that 34% 
of the mass in the Universe is concentrated in filaments, as the 
numerical results of [Hoffman et all (12012l) indicates. We then 
have that 


are clearly different: in particular, for a typical galaxy cluster 
with L c i - 1 Mpc at z\ = 0.5 and for a source at g s = 1 we find 
that the coupling constant between the shear y and the tidal field 
is about -20C. As a result, in this context (for a galaxy cluster) 
we expect that the intrinsic alignment produced by the instanta¬ 
neous tidal field is negligible. 

3.2. Particular cases 

We now consider an external Keplerian field, that is, the poten¬ 
tial generated by a second galaxy of mass M ext , assumed to be 
spherical and centered at R. The tidal acceleration and the tidal 
potential acting on a star of the first galaxy are 


A(M C j) ^ 0.34p n ,4 = 4.67 x lO 7 M 0 pc _I . (45) 

If we place the filament in such a way that it crosses the source 
plane at R = (R, 0), with R = 500 kpc, and use Eq. (l34l >. we 
obtain (Fig. [2} 

e s = +———=; +1.4x 1(T 3 . (46) 

5cr 2 R 2 T\ - Ti 

We observe that in this case the dependence on the distance is 
e s = O (R- 2 ). 


4. Discussion 


&i, tidal, M — 
^Aidal,M = 


GM ext 

R 3 

GM ext 

2/7’ 


' R i R i\ 

(37) 

' R i R i\ 

Av-3 y*i. 

(38) 


If we place the external galaxy at R = (R. 0,0), using Eq. <im 
we obtain (Fig.0 


s= 3GM ext n 
+ 20cr 2 R 3 Ti-T 2 


+4.6 x 1(T 6 , 


(39) 


where in the last step we used M ext = 10 11 M e and R = 500 kpc. 
We note that s s = 0 (R 3 ). 

We now examine a different case: the external field is gener¬ 
ated by a dark matter (DM) filament, and the galaxy is outside it. 
We are considering a tidal field generated by an (infinite) linear 
distribution of matter along the line of sight X 3 


&i, tidal, /l 




^tidal,/l — + 


GA 

R 2 


‘ R 2 




(40) 

(41) 


4.1. Comparison with literature 

It is common (e.g., iJoachimi & Bridlj 120101: iJoachimi et al.l 
1201 ll) to consider the galaxy number density contrast-intrinsic 
ellipticity correlation, that is, 

Cgi(x) = <(5„(x 0 )s+(x 0 + x)) Ao , (47) 

where the correlation function only depends on the modulus of 
x because isotropy, 6 „ is the galaxy number density contrast, and 
is the real part of the complex (intrinsic) ellipticity in a co¬ 
ordinate system with the real axis along the line connecting the 
two correlating points Xq and x, as in Eqs. (|5}, (( 6 )i, and ([TJ. In 
contrast to Eqs. ([5} and following, the correlation in Eq. ([47} is 
carried on in the three-dimensional real space. We now define 
the Fourier transform / (k) of a function /(x) as 

/(k) - J /(x)e' k x d 3 x . (48) 

4 An important point to keep in mind is that the large-scale structure is 
strongly time-dependent and that the mean high-mass-cluster distance 
changes with time dSpringel et alJl2005l) . Here, we are considering the 
local Universe, that is, regions where the redshift is much lower than 
unity. 
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Since 

<Wk) = b e 6 (k) , (49) 

k 2 -k 2 

el(k) = - C(k 2 - k 2 )U( k) = 4nCG^ r ^p cr Qj(k) , (50) 

where b s is the galaxy bias and in the second equality of Eq. (l50l > 
we used Eq. dm the Poisson equation and the relation p(k) = 
p cr Q ln h(k) (we neglect the three-dimensional Dirac delta in the 
origin). If we introduce the elliptical galaxy fraction f e \\ (because 
only elliptical galaxies align in our model), we find 

(2/r)V 3) (k - k')P gl (k,z) = { 6 n (k)sl(k')) 

- (2;r) 3 d (3, (k - k')4nb & f e nCG Pcl klMk, z) , (51) 

where P^(k,z) is the density contrast power spectrum at the 
galaxy redshift z and where we have approximated the term 
(k 2 - kz; )/k 2 to unity. We may now consider Eqs. (6) and (19) 
of lJoachimi et al. ( 201 11. that is, 

P< “ AB A«= -^*W) PAk - l) (rrsj (£f ■■ (52) 

where D(z) is the linear growth factor normalized to unity today, 
//> is the (nonlinear) power spectrum of the density contrast, z 
and L are the redshift and the absolute luminosity of the early- 
type galaxy, and the arbitrary reference values are z.n = 0.3 and 
Lo, the latter corresponding to an absolute r-band magnitude of 
-22, passively evolved to z — 0. The parameters AC i, //, and 
/j have been measured by iJoachimi et al.l ( 2011 1 for comoving 
transverse separations greater than 6/z 1 Mpc, 

AC iP cr = 0.077 ± 0.008 , (53) 

yS= 1.13!“;“, (54) 

V = -0.27!°;*° . (55) 

Comparing Eq. ( |5H to Eq. (l52l > and neglecting the redshift and 
luminosity dependence, we obtain that Eq. (l53l > corresponds to 
the term 

-47rCGp cr ^ 3.2 x 1(T 6 (56) 

of Eq. (ITTl i. where we set / e n = 1 because lioachimi et alJ (1201 ll) 
considered only early-type galaxies for the shape measurements. 
The result is clearly inconsistent. It is unlikely that this inconsis¬ 
tency is due to our use of the simple King model because the dif¬ 
ferent approach adopted in Appendix |Bj where we do not make 
use of the stellar distribution function, yields similar results (cf. 
Eqs. <[39]i. ©. C12Di and (IB. 101 0. Then, at least at large scales 
( r p > 6/z 1 Mpc), our model does not reproduce the observations. 

The incompatibility of Eq. (l56l) with Eq. ( 1531 ) means that 
the galaxy deformation due to the instantaneous external tidal 
field cannot yield the observed IA signal. A possible explana¬ 
tion is that the IA signal is caused by the galaxy formation pro¬ 
cess and/or its merging history. To obtain analytic results, these 
processes therefore need to be linked to the external tidal field. 
In particular, at least for the merging history, the velocity shear 
field needs to be con s idered because it was recently discovered 
dHoffman et alJl20l2t iLibeskind et al.ll2014t iLee & Choil 120 1 4l> 
that mergers preferentially occur along the velocity shear minor 
eigenvector. A detailed analysis of this process is beyond the 
aims of this paper. 


4.2. Halo case 

Until now, we have only considered the galaxy with its total (lu¬ 
minous plus dark) mass. We have stopped the integration of the 
luminous quadrupole at r lnax < r,, . However, the ellipticity of 
the whole DM halo remains unclear. We remark that it is not 
obvious that Eq. (0 is the correct choice for the halo elliptic¬ 
ity because the halo luminous quadrupole cannot be measured in 
real surveys. For the same reason, it is unclear how to choose 
r max , but for the condition |f/ 0 (r max )| » |f/ t id a i(r max )|. Still, it is 
interesting to consider the extreme case of a dark matter halo 
M — 2 x 1O 13 M 0 , r a = 560 kpc and cr„ = 250 km s -1 , as found 
bv lGavazzi et al.l d2007h from an ensemble average over 22 mas¬ 
sive halos. Even now, we obtain (taking r corc = 0.01 ■ r u and 

A nax = 0.80 ■ r tl ) 

Ch a io - -1.3 x 10 17 yr 2 , (57) 

-47rC ha i 0 Gp cr = 0.0042 , (58) 

which is still almost a factor ~ 20 lower than the value of 
Eq. (l53l >. that is, even if instantaneous tidal fields produce a 
greater IA on dark matter halos than on galaxies, the effect is still 
much too weak to account for the observed signal. Arguably, the 
galaxies observed bv lToachimi et al.l ( 2011 1 are, on average, less 
massive objects, and therefore the value of Eq. ( l57l > should be 
taken as an upper bound. In addition, we stopped the integra¬ 
tion near the halo boundaries (r max = 0.80 ■ r tl ), and the resulting 
ellipticity is therefore greater than that we would obtain with a 
smaller r max . This shows that our model is essentially inadequate 
to explain the observed amount of IA. 


5. Conclusions 

The exact origin of intrinsic alignment of galaxies is still unclear, 
and there are no analytic estimates of the amount of IA caused 
by different possible processes. We here support the standard 
view on IA (i.e., it is caused by formation and accretion pro¬ 
cesses) by a reduction ad absurdum. To arrive at this result, we 
estimated the amount of IA in elliptical galaxies due to an ex¬ 
ternal instantaneous tidal field by considering the galaxy stellar 
distribution function and the luminous quadrupole. In addition, 
in Appendi x [A] we determined the typical oscillation time-scale 
for an elliptical galaxy modeled as a rigid body, subject to an ex¬ 
ternal tidal field, and in Appendix[B]we determined the ellipticity 
of a galaxy (subject to an external instantaneous tidal field) in a 
different way by studying its equipotential surfaces. The main 
results are the following: 

1. The distribution function approach allows us to analytically 
determine the dependence of the intrinsic ellipticity on the 
tidal field, heuristically formulated by ICatelan et al.l (1200 lb . 
in terms of the properties of the galaxy (i.e., its mass, 
size, velocity dispersion, and stellar distribution function), 
Eqs. (l34l > and (l35l >. 

2. The intrinsic alignment signal obtained when our model 
is applied to an elliptical galaxy is negligible with respect 
to the observed one (at least at large scales, > 6/z 1 Mpc, 
IJoachimi et al.ll2()TTl) . cf. Eqs. (l53l) and d56t . Thus one has 
to consider the galaxy formation process and/or its merging 
history. 

3. When our model is applied to the whole galaxy halo, the 
intrinsic alignment signal increases (but it is still inconsistent 
with the observed signal), Eq. d58l . 
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The work we described here is a step toward a simple phys¬ 
ical understanding of the bias introduced in weak-lensing mea¬ 
surements by IA. In the past, IA has been regarded as a source 
of concern for cosmic-shear measurements. More recently, it has 
been regarded as an opportunity to investigate the physical prop¬ 
erties of galaxies, their DM halos, and their formation history. In 
this perspective, it is important to develop analytic models of IA, 
such as those presented here, to be able to interpret the results of 
future weak-lensing surveys. In this respect, a complete model 
of IA could even be used directly to reconstruct the local tidal 
field acting on elliptical galaxies. 

The theoretical analysis could be improved by applying the 
techniques we presented (i.e., the stellar distribution function 
method) to the study of the formation and accretion processes, 
in order to determine their contributes to the IA signal. It would 
also be interesting to study the contribution to IA due to contin¬ 
uously applied tidal fields (instead of only considering the tidal 
field acting at the emission epoch, as we have did here), and the 
evolution of IA with redshift. 
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Appendix A: Rigid-body approximation 

In this appendix we compare the time scale of motions of galax¬ 
ies modeled as rigid bodies with the time scale of their internal 
dynamics. In particular, we consider the (rigid) oscillations of el¬ 
liptical galaxies and the precession of spiral galaxies. We show 
below that internal dynamics is much faster than rigid-body dy¬ 
namics, thus confirming the results of ICiotti & Duttal (1 994 ) for 
elliptical galaxies. 


Appendix A. 1: General case 

We assume that the distribution of mass of the galaxy only de¬ 
pends on the elliptical radius r and that its principal axes are 
aligned with the coordinate axes: 

p(x, x, t) = p(x) — p(r), (A.l) 


x\ 


r\ 


X 2 


— + — + — 


ri 


x-i 


r 3 


(A.2) 


where r\, r 2 and r 3 are the semi-axes of the ellipsoid (p - 0 for 
r > 1). Then the inertia tensor is 


I U = 


= J' P(x)(||x||% - XiXj) d\ 


(A.3) 


where 6 is the Kronecker delta, and we placed the center of mass 
of the galaxy in the origin. With this definition, the moment of 
inertia /ft along the unit vector n is 


Ifi = h l \ijh ] . 


(A.4) 


For an elliptical mass distribution whose principal axes are par¬ 
allel to the coordinate axes we obtain 


I = Tor x r 2 r 3 diag (r? + rj, r\ + rj, r\ + r 2 ) , 


^ 4-r f 1 

n -~J„ 


p(r)r 4 dr . 


(A.5) 
(A. 6 ) 


The potential energy W of a body subject to the action of the 
external field U ext (x) (whos e La placian is null in the region oc¬ 
cupied by the galaxy) is dJacksonlll998l) 


/ 


W = p(x)t/ ext (x)dA 


MU ext ( 0 ) + -Q.jffd’UUO) 


Qu = 


Jp(x)(3x,xj - ||x|| 2 6 !7 ) d 3 . 


(A.7) 
(A. 8 ) 
(A.9) 


where M and Q t j are the mass and the quadrupole tensor of the 
body. For an elliptical mass distribution with principal axes par¬ 
allel to the coordinate axes we obtain 


<3 = Tor\r 2 r 3 diag ( 2 r 2 - r\ - r\, 

- r 2 + 2r\ - r 2 , -r 2 - r 2 + 2 rfj . (A. 10) 
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Fig. A.l. Spherical coordinate system adopted in the text, Keplerian 
(M ext ) and DM filament cases. 


Appendix A.2: Elliptical galaxies: oscillation period 

For a Keplerian potential, the galaxy energy is 


Wm = ~ 


GM ext M 

R 


GMmQij I R,Rj \ 

""oM S 5 " T 


(A.l 1) 


Without loss of generality, let fi be diagonal and let 6 and <p be 
the angles between the x 3 axis and R and between the xj axis and 
R-R-X3, in a spherical coordinate system (see Fig. lA. It . In other 
words, we consider a coordinate system whose axes are parallel 
to the principal axes of the galaxy. Having kept fi diagonal, we 
have to rotate the tidal tensor. Considering that Tr fi = 0, we find 


Wm = 


GM„,M GM. 


ext 


(fin cos 2 * (^) sin 2 ( 0 ) 


R 2 Ri 

+ Q22 sin 2 (<p) sin 2 ( 0 ) + <233 cos 2 ( 0 )) . (A. 12 ) 


We can in principle determine the oscillation equations of the 
galaxy via the Euler-Lagrange equations, using the correct ki¬ 
netic energy expression. This approach leads to complex equa¬ 
tions and therefore we assume that (i) the elliptical distribu¬ 
tion of mass is prolate (i.e., r\ — r 2 — r eq < r 3 ), from which 
hi = 1 22 = leq and fin = fi22 = fieq, and (ii) there are no initial 
intrinsic rotation (<p(0) = 0 ) or precession motions; there could 
be only nutation (6(t) 4 0). Using the Euler-Lagrange equations, 
we then obtain 


tude, we find that the period t 0sc ,m of an oscillation is 


. R 3 

lose ,M — 2 71 -» I 

V GM ext 
=* 3.0 x 10 11 x 



T/ e\t 

10 n M o 


yr. 


(A. 14) 


We remark that Eq. (IA. 1 4b is reminiscent of the expression of 
the free-fall time oc (Gp )~'/ 2 oc f oscM x (D/2R) 3 ! 2 , being D 
the diameter of the galaxy. Using the values adopted here for the 
external mass M ext = 10 11 M G and for the distance R = 500 kpc, 
the period is f osc m - 1.0 x 10 11 yr, which is greater than the age 
of the Universal (the numerical integration of Eq. (1A. 1 3b is in 
accordance with the approximate result). 

This period has to be compared with the time scale on which 
the shape of a galaxy changes, because of the external gravita¬ 
tional field. A rough estimate of this is the time necessary for a 
star to travel across the galaxy (e.g jFleck & Kuhnll2003l) 


Across,ell — 


D 


(A.15) 


where cr„ is the stellar velocity dispersion and D is the diame¬ 
ter of the galaxy. For D - 40 kpc and cr„ - 200 km s we 
obtain f C ross,eii — 2 x 10 8 yr <s t 0 sc ,m- In other words, an ellipti¬ 
cal galaxy, not aligned with the gravitational tidal field, deforms 
itself before completing a rigid oscillation. 

If we consider a DM filament along the xj axis, using 
Eq. (ED, we can write the galaxy energy as 


W A = +2MGA log 



GAQ U 
3 R 2 




(A.16) 


where £' is an arbitrary constant introduced to have an adimen- 
sional logarithmic argument, and 6 = diag(0,1,1). As before, 
we rotate the coordinate frame in such a way that fi,y be diag¬ 
onal (see Fig. IA. 1 1 . we use spherical coordinates and assume 
that the galaxy is prolate, that there are no proper rotation nor 
precession motions, and that the position of the galaxy is fixed. 
We then obtain 

GA 

M = J^2 sin( ’ 26 ') • ( 3 eq - Q33) ■ (A. 17) 

Using the same approximations as before and Eq. (l45l >. we may 
estimate the period of a small rigid oscillation of the galaxy 


M sin ( 20 ) ■ (^q - 3 m) . (A. 13) 

We explicitly note that this equation is not self-consistent: in 
fact, in principle it should be completed with the equations of 
the relative motion between the galaxy center of mass and the 
monopole. However, we are not interested in a complete treat¬ 
ment of the problem, but only in an analytic estimate of the time 
scale of an oscillation. Therefore we assume, without clarifying 
the physical mechanism that would permit this, that the relative 
position between the external monopole and the galaxy is fixed, 
and that only the orientation of the galaxy can vary. With these 
assumptions, Eq. (IA. 13b is sufficient to describe the dynamics of 
the system. To obtain the time scale of an oscillation, we make 
the small-angle approximation (0 « 1 ) and recall that <333 > fi eq 
(prolate galaxy), obtaining the pendulum equation. Considering 
that the components of I and fi are of the same order of magni- 


t° sc ,A - 2nR y 2G ^ 

^ 1.0xl0 10 x(—lyr, (A.18) 

\ 1 Mpc / 

which in the case considered throughout the paper (R = 500 kpc) 
corresponds to 5.0 x 10 9 yr. This is shorter than the period of 

5 It is interesting to consider the case of stellar globular clusters: as¬ 
suming a distance from the host galaxy of R = 30 kpc, we obtain a 
period t osc>c i = 1.5 x 10 9 yr. We emphasize that this value has to be 
increased by the multiplicative factor -^/l eq /(^33 ~ fieq)- 

6 This procedure is more complicated than in the Keplerian case be¬ 

cause d is not a three-dimensional delta. It can be written as <5 ,7 = 6ij - 

ic it 1 Xj_ 1 , where X| = (1,0,0) is the versor of the first axis, parallel to the 

DM filament. When we change the coordinate frame, we have to rotate 

R and Xi; we then obtain W A (d,<f>) = const - 3 -1 G/l/?~ 2 (Qn cos 2 0(1 + 

sin 2 6) + Q 22 sin 2 0(1 + sin 2 6) + <3 33 (1 + cos 2 8)) . 
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oscillation found in the Keplerian case, but nevertheless greater 
than f C ross,eii- We therefore have to drop the rigid-body approx¬ 
imation and deepen the description of elliptical galaxies to ac¬ 
count for the internal degrees of freedom, as done in Sect. [3] and 
Appendix |B1 


Appendix A.3: Spiral galaxies: precession period 


The key feature that distinguishes spiral galaxies from ellipti¬ 
cal ones is the dominance of ordered motions over chaotic ones, 
that is, the characterizing presence of the angular momentum. If 
we place a rigid body with an angular momentum in an exter¬ 
nal field, it starts precessing. In this subsection we estimate the 
precession period for a spiral galaxy, taken as a rigid body. 

The precession period f prec of a rotating rigid body is 


t 


prec — 


27 rLsin( 0 ) 
T (9) 


(A. 19) 


where L is the angular momentum of the spiral galaxy, 6 is the 
angle between L and the external force, and t( 6 ) is the amount of 
the momentum of the external force. In spiral galaxies the lumi¬ 
nous mass distribution is very different from that of the dark mat¬ 
ter ; in our calculation we can detach the two contributions be¬ 
cause the DM distribution is (approximately) spherical and does 
not generate a torque on the visible mass, in which we are in¬ 
terested. To determine the total stellar angular momentum L*, 
we can use the rigid-body formula L* = 1 * 3307 *. In spiral galax¬ 
ies stars at different distances from the center of the galaxy have 
different angular velocities (to* = to*//), where ( is the distance 
from the symmetry axis of the spiral). Therefore it is sensible to 
use a weighted angular velocity to*, defined as 


W* - L*/l*33 


-I 


/4*(/)I*(4d/ I fSJ£)d£, 




(A.20) 


where Z*(f) is the projected stellar surface density at distance 
£ from the spiral symmetry axis and o*(T) is the stellar tan¬ 
gential velocity. In spiral galaxies, the stellar tangential veloc¬ 
ity is approximately constant, v+it) - 0 *, and the stellar col¬ 
umn density follows an exponential law with length scale 4 > 
£*(4 = So expf-fyfo). We then obtain 


_ v+ Jg 00 1 2 exp(-l/4) d£ Vir 
/“^ex p(-£/£ 0 )d£ 34' 

Using this expression in Eqs. (IA. 13b and ( IA. 17k we find for an 
oblate galaxy (r 3 < r\ = r 2 = r eq ) 

_ 27ro* I *33 R 3 1 

pre c.M - 3G4 __ Mex( cos(0) > 

_ 7TP.fr 1*33 R 2 1 
prec - 1 " G4<2*eq-0*33 A COS(d) ' 


(A.22) 
(A.23) 


Assuming p* = 200 km s 1 and 4 = lOkpc, we obtain 

fprec ,m = 1-2 x 10 13 ^ (cos(6>)) yr , (A.24) 

fprec, /t = 2.7 x 10 10 —!%-(cos^f 1 yr . (A.25) 

Therefore the time scale is longer than or similar to the age of 
the Universe, and it is also longer than the deformation time of 
spiral galaxies (similarly to that of the elliptical galaxy). 


Appendix B: Equipotential approximation 

In Sect. [3 we have obtained the expression of the intrinsic el- 
lipticity of an early-type galaxy subjected to an external tidal 
field. To do so, we have calculated the luminous quadrupole, 
making use of the stellar distribution function. In this appendix 
we present another approach, which is less complete but has the 
advantage of having a clear and simple physical understanding. 
In particular, we model the deformation of the galaxy by means 
of the equipotential surfaces of the total gravitational potential. 
In this approach, we “assume[s] that the local galaxy density 
is produced approximately by stars near their zero-velocity sur¬ 
faces” (ICiotti & Duttall 19941) . As in Sect. [3 (i) we start with an 
unperturbed spherical galaxy, (ii) we take an external tidal poten¬ 
tial and add it to the (unperturbed) galaxy potential, thus ignor¬ 
ing the changes in the gala xy potential induced by the deforma¬ 
tion o f the galaxy (see also ICiotti & Duttall 1994tfBertin & Vairil 
2008) , and (iii) we assume that the galaxy immediateyl reacts 
to a change of the external gravity field by modifying its shape 
accordingly (see Sect.[3for details). 


Appendix B. 1: General case 

In the absence of external fields, the galaxy potential Uq obeys 
the Poisson equation 

V 2 f/ 0 = 47rGp 0 , (B.l) 

where po is the unperturbed galaxy mass distribution. We now 
introduce a (weak) external potential U ext , so that Uq 1 —> U - 
Uq + I/ext- The introduction of the external field changes the 
equipotential surfaces of the galaxy. Given a volume 'Vo en¬ 
closed by a particular equipotential surface <YVq - {Vx|I/o(x) = 
/-’(,} at energy Eq of the unperturbed potential, we consider the 
corresponding equipotential surface d' V for U 

£fV = {Vx| U(x) = t/ 0 (x) + U ex,(x) = £0 + 6E = E) . (B.2) 


The energy shift 5E is chosen in such a way that the mass inside 
the surface does not change: 


J' &(Eq - t/ 0 (x))p 0 (x) d 3 r = J' ©(£ - U(x))p 0 (x)d 3 x . 


(B.3) 


A Taylor expansion of the right-hand side of Eq. (IB .3b for low 
6 E - U e xt gives 


J 6{Eq - Uo(x))p 0 (x){ 6 E - f/ ext (x)) d 3 x = 0 . 


(B.4) 


Since p 0 has spherical symmetry, po and Vf/o are uniform on 
(YVq , and we finally obtain 


6 E = 


£y„ t/ext(x)po(x)|Vf/ 0 (x)| 'd 2 r £y„U ex t(x)d^t 


Wo 


&V Po(x)|Vf/ 0 (x)p'd 2 x 


iv„ d2 *' 


(B.5) 


If {/ext = Utidai = 4 'ijXiXj, that is, the external potential is a tidal 
one, 6 E is equal to zero, because <I/ ; is traceless and <YVq has 
spherical symmetry. 

We use the galaxy logarithmic potential of Eq. (l28l >. If we 
introduce an external tidal potential, the equipotential surface 
becomes an ellipsoid (see Fig. ED- We place the galaxy center 
at the origin and align the coordinate axes along the eigenvec¬ 
tors of the tidal tensor (h, 7 . Then, we can evaluate the deviation 
from the circular shape of a particular equipotential surface with 
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Fig. B.l. Schematic representation of the deformation of a two- 
dimensional contour of an equipotential surface. The dotted line is the 
unperturbed contour, the solid line the perturbed contour. In the tidal 
approximation, the variations along the positive and negative directions 
are the same. 


The values obtained through the equipotential approxima¬ 
tion (IB. 81 ) and (IB. 101 ) are higher than those obtained through 
the distribution function method ( l39l ) and (l46l) : the reason is 
that with the equipotential approximation we only consider an 
outer isopotential surface, which is more deformed than the in¬ 
ner ones. Instead with the distribution function method, more 
realistically, we are “weighting the ellipticities of the isophotes” 
of the galaxy from its center to r max . Equations (IB. 8 b and (IB. 10b 
have the same dependence on R as was found with the distribu¬ 
tion function method (see Fig. E](. 


radius r max < r tr by expanding to first order its semi-axis varia¬ 
tions 6 j in Eq. (|B3 . 


6i = 


"LTmaxTr 


MG(r a + 

Lore)/(Lnax + Lore) + ^^max^tr^ 


(B. 6 ) 


We kept for generality 6 E, even if it vanishes for an external tidal 
field. Equation (IB. 6 b allows us to compute the intrinsic ellipticity 
of a particular isophotal of a galaxy subject to a tidal field as 
observed along any direction. For this, we just have to project an 
ellipsoid with semi-axes r max + 5i, r max + 62 and r max + (>3 along 
the line of sight. For example, if the ellipsoid is observed along 
X3, we would have from Eq. (0 


Or 

max 


(B.7) 


Appendix B.2: Particular cases 

For the Keplerian tidal field ( 138b . placing the external spherical 
galaxy at R = (R, 0,0), using (IB. 7b . and assuming r max , r tr <sc R, 
we obtain (Fig. Ell 


3 M ex t T'max^tr ^ max + f core 

4 Ad R ^ Hr + y core 


1.6 X 10 5 , 


(B. 8 ) 


where in the last step we used M ext = 10 11 M Q , r max = 10 kpc and 
R = 500 kpc. We note that s s = 0 (R 2 ), and that the ellipticity 
becomes lower for inner equipotential surfaces. 

For an external DM filament directed along the line of sight 
and distant R = 500 kpc from the galaxy, if we assumtfl 


M—— » 2 ? (B.9) 

f max F core ^ 

we obtain from Eqs. (TiTl i. (l45l > and (IB. 7b (Fig.E} 


'Ir max r tr r, nax + r cole 
MR- /'n + r core 


^4.9x 10 ^ 3 . 


(B.10) 


We observe that in this case the dependence on the distance is 
e s = 0(R 2 ); again, inner equipotential surfaces have lower el¬ 
lipticity. 


7 With the values used in this paper, 2/b' max r t 2 r (r core + r max ) ^ 
10 - 2 MR 2 (r COK + r tr ), and our assumption is justified. 
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